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ABSTRACT 

m 

O 

^ . X-ray spectra of the late-type star AB Dor, obtained with the XMM-Newton 

■ satellite are analysed. AB Dor was particularly active during the observations. 

. An emission measure reconstruction technique is employed to analyse flare and 

quiescent spectra, with emphasis on the Fe XVII 15 — 17 A wavelength region. 
The Fe XVII 16.78 A/15.01 A line ratio increases significantly in the hotter flare 
plasma. This change in the ratio is opposite to the theoretical predictions and is 
attributed to the scattering of 15.01 A line photons from the line-of-sight. The 
q , escape probability technique indicates an optical depth of ~ 0.4 for the 15.01 A 

line. During the flare, the electron density is4.4+ 2 ;£xl0 10 cm and the fractional 
! Fe abundance is 0.5 ±0.1 of the solar photospheric value (Grevesse et al. 1992). 
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Using these parameters a path length of ~ 8,000 km is derived. There is no 
evidence for opacity in the quiescent X-ray spectrum of the star. 
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1. Introduction 

Opacity in the solar corona has often been invoked in order to explain the discrepancies 
in the observed and theoretical spectra of several spectral lines, including O VIII, Fe XVII, 
Ne IX and Mg XI (Acton 1978; Schmelz et al. 1997). During the process the photon of a 
resonance line is absorbed by an ion in the ground state, and this excited ion emits a photon 
at the same wavelength but in a new direction. This process of absorption and re-emission 
is called resonance scattering. The optical depth at line center can be written as (Mitchell & 
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Zemansky 1961): t oc \fij\J y^^^T 71 ^ wnere A is the wavelength in A, the oscillator 
strength of the transition, M the mass of the absorbing atom in atomic mass units, T the 
temperature in K, / the path length in cm and n- u n e \, riu, n c the number densities (in 
cm -3 ) of ions in the lower level i, element, hydrogen and free electrons, respectively. Some 
emission lines can have sufficiently large oscillator strengths to lead to high optical depths. 
Depending on the density and geometry of the emitting region, photons can be scattered 
out or into the line-of-sight. The method used to study coronal opacity (originally applied 
to chromospheric line of C I; Jordan (1967)) relies on pairs of lines from the same element 
and ionization stage but with different oscillator strengths. The line with the low oscillator 
strength, which is used as a reference, is assumed to be optically thin and is often a forbidden 
transition. The analysis of solar active regions by Schmelz et al. (1997) revealed statistically 
significant opacity only in the Fe XVII line at 15.01 A. 

Based on the EUVE observations of cool stars, Schrijver et al. (1994) suggested that 
resonance scattering could account for the low line-to-continuum ratio seen in their coronal 
spectra. This suggestion was later ruled out by Schmitt et al. (1996) who demonstrated that 
the count fluctuations in the EUVE spectrum suggest that the bulk flux arises from line 
and therefore the continuum emission had been significantly overestimated (see also Jordan 
(1996)). Laboratory experiments carried out with the LLNL EBIT-II, have indeed shown a 
large number of newly identified transitions in the 60 — 140 A wavelength range that add up 
to a significant flux (Lepson et al. 2002). 

Analysis of HST and FUSE observations of cool stars have shown significant deviations 
from the optically thin approximation in the transition region lines of C III, C IV, Si IV and 
O VI (Mathioudakis et al. 1999; Bloomfield et al. 2002). The wealth of high quality X-ray 
observations from the Chandra and XMM-Newton missions has allowed the application of line 
ratio techniques to the coronae of cool stars. These techniques were applied to the quiescent 
coronae of the RS CVn binaries Capella (Phillips et al. 2001) and II Peg (Huenemoerder 
et al. 2001) and revealed no evidence for appreciable opacity. A systematic investigation 
of opacity effects in 26 late-type stars led to the same conclusion (Ness et al. 2003). More 
recently Testa et al. (2004) claim the detection of opacity in the O VIII lines of II Peg. It 
therefore seems that, while there is strong evidence for resonance scattering when observing 
individual solar active regions and flares, these effects are not very well pronounced when 
observing the emission integrated over the stellar surface. 

AB Dor (HD 36705) is a young (20 — 30 Myr), active and rapidly rotating (vsini = 
90 km s _1 ) K0 dwarf at a distance of 14.9 pc (Unruh et al. 1995). In the present paper 
we analyse spectroscopic X-ray observations of AB Dor obtained with XMM-Newton. Our 
particular interest is on the study of resonance scattering in a flare by means of Fe XVII line 
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ratios. 



2. Observations and Data Reduction 

The observations presented in this work have been obtained with XMM-Newton's Re- 
flecting Grating Spectrometers RGS1 and RGS2 (den Herder et al. 2001). Data were reduced 
with the XMM Science Analysis System (SAS) software v5.3.3, updated with the latest cal- 
ibration files. The light curve derived from the EPIC pn data is shown in Fig. 1, where 
we can see that two flares were observed. Using the EPIC pn light curve as a guide, we 
have extracted RGS spectra for the first flare and quiescent time intervals (FLARE and 
QUIES1+QUIES2, respectively) as indicated in Fig. 1. Based on the analysis of the back- 
ground light curve extracted from CCD9 of RGS2, time intervals which show high count 
rates, probably due to proton events, were excluded. The RGS1 and RGS2 spectra are 
divided into 3400 channels with widths ranging from 0.007 A to 0.014 A and cannot be 
summed directly as their wavelenghts grids do not coincide. Spectra were re-binned into a 
new grid with a constant width of 0.02 A before the analysis. Consequently, we have defined 
a total Line Spread Function (LSF) for the summed spectrum, which is the sum of the single 
LSFs weighted by the effective area (see Scelsi et al. (2004) for more details). 

The analysis performed on the RGS spectra is based on the plasma emission measure 
distribution (EMD) versus temperature reconstruction technique. We have used the IDL (In- 
teractive Data Language) based software Package for Interactive Analysis of Line Emission, 
PINTofALE vl.5 (Kashyap & Drake 2000) for this purpose. The analysis can be summarized 
in the following steps: 

identification and fit of the lines. For each spectrum we have identified the strongest 
lines using the CHIANTI v4.0 atomic database (Dere et al. 1997; Young et al. 2003) together 
with the ionization balance of Mazzotta et al. (1998); as there is no analytical expression 
for the total LSF, we have approximated the profile of each line with a Lorentzian having 
a width fixed at the value which provides the best match with the actual total LSF. Scelsi 
et al. (2004) have demonstrated that this is a good approximation and that most of the 
flux measurements performed in the total spectrum are accurate. It is worth noting that 
since the wings of the line profile are extended, in the spectrum there is the formation of a 
"pseudo-continuum" which makes measurements of the line flux difficult; multicomponent 
fitting involving the adjacent transitions was therefore necessary for most of the lines; 

EMD and predicted spectrum. We have reconstructed the emission measure distribution 
from measured line fluxes using the Markov-Chain Monte Carlo (MCMC) method described 
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by Kashyap & Drake (1998). In order to derive an EMD we have chosen intense lines, i.e. 
those for which the Lorentzian was a good approximation of the actual LSF and did not 
fall on the bad pixels of the detector. In the lines selected the Lorentzian profile should 
reproduce at least 95% of the actual flux. The chosen lines are reported in Table 1. Finally, 
using the EMD we obtained a predicted spectrum. 

3. Results and Discussion 

In Figures 2 and 3 we show the observed and predicted spectra obtained from the 
analysis of the flare and quiescent data, respectively. While the observed over predicted 
ratio for the Fe XVII 16.78 A line remains the same in the flare and quiescent spectra of AB 
Dor, the corresponding ratio for Fe XVII 15.01 A is reduced significantly in the flare. 

There are some known atomic data problems regarding the Fe XVII 3d — > 2p and 
3s — > 2p lines located at approximately 15 and 17 A. The observed intensity ratios of the 
3s — > 2p lines (17 A) to the 3d — > 2p lines (15 A) are consistently larger than those calculated 
by collisional-radiative models. There have been some experimental attempts to explain this 
discrepancy. Laming et al. (2000) measured the 3s/3<i ratio using the NIST EBIT and found 
that their data are in agreement with the ratios predicted using a distorted wave collisional 
model. A similar experiment was carried out by Beiersdorfer et al (2002) at the LLNL 
EBIT-II: the comparison between their data and theoretical calculations shows that resonant 
excitation makes a significant contribution to the 3s — > 2p lines intensities. This suggests 
that additional processes have to be taken into account in the interpretation of astrophysical 
plasmas. Doron & Behar (2002) and Gu (2003) have shown that the contribution of radiative 
recombination (RR), dielectronic recombination (DR), resonant excitation (RE) and inner- 
shell collisional ionization (CI), in addition to collisional excitation (CE) can rectify some of 
the disagreements between theory and astrophysical observations. Both studies show that 
as the temperature increases the 16.78 A/15.01 A ratio decreases, since the net relative 
contribution of DR, RR, RE and CI processes for the 3s level, compared with collisional 
excitation, decreases with increasing temperature. These new atomic data of Fe XVII have 
been used in our analysis (Doron & Behar 2002). 

The plasma becomes hotter during the flare. The EMD becomes steeper (EMD q oc 
T L5±0 - 7 , EMD/ oc T 2 - 7±L0 in the temperature region log T = 6.5-7.1) with the minimum 
shifted to higher temperatures (from log T q = 6.5 to log Tf = 6.7). We have folded the 
Fe XVII 15.01 A and 16.78 A line emissivities (log T pea k = 6.75) with the quiescent and flare 
EMDs. The contribution function peaks at log T q = 6.8 and log Tf = 6.9 for the quiescent 
state and flare respectively. The Fe XVIII 16.073 A to Fe XVII 16.78 A ratio is often used as 
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a temperature diagnostic, (Phillips et al. 1997), and applied to our observations gives log T q 
= 6.6 and log Tf = 6.8. Finally, a fit of the EPIC pn observations of the same flare carried 
out by Giidel et al. (2001) requires a higher temperature compared to the quiescent state. 
Although the temperature is therefore higher, the observed 16.78 A/15.01 A ratio increases. 

In Table 2 we list the observed Fe XVII 15.01 A and 16.78 A line ratios for the quiescent 
state and flare of AB Dor. The theoretical ratios of Doron & Behar (2002) and CHIANTI v4.0 
are also listed. While for the quiescent observation the ratio is consistent with the theoretical 
values, as the temperature increases during the flare, the ratio increases significantly: this 
observed behavior is opposite to the theoretical one. Using the APED line list (Smith et al. 
2001) we have folded the emissivities through the emission measure distributions and LSF 
and found that line blending could account for up to 7% and 2% of the 15.01 A and 16.78 A 
line flux respectively. We attribute the increment of the Fe XVII 16.78 A/15.01 A line ratio 
that is seen during the flare to the resonance scattering of Fe XVII 15.01 A photons out of 
the line-of-sight. 

Using the line ratio method which has been applied successfully to both stellar and solar 
spectra (Mathioudakis et al. 1999; Bloomfield et al. 2002; Schmelz et al. 1997), the reduction 
in the 15.01 A line flux allows us to derive the path length of the scattering layer. The photon 
escape probabilities as a function of optical depth are described in Kastner & Kastner (1990). 
Here we assume an inhomogeneous geometry where the emitters and absorbers are spatially 
distinct. The ratio between the fluxes of the two emission lines can be written as: 

-^16. 78 _ -Pl6.78( r o) £l6.78 ^ 
-^15.01 -Pl5.0l( T b) £l5.01 

where -F\ 6 .7 8 and Fi 5 . i are the observed fluxes, -Pi6.78( r o) and Pi5.oi(to) are the pho- 
ton escape probabilities and £i6.78 and £15.01 are the emissivities for the optically thin 
(16.78 A) and optically thick (15.01 A) line respectively. The ratio of the observed fluxes is 
-F16.78/-F15.01 — 0.71, while the ratio between the emissivities is £i6.7s/£i5.oi — 0.54 (Doron 
& Behar 2002). If we assume that -Pi6.78(t"o) — 1, from eq.(l) we have -Pi5.oi(fo) = 0.76 
which gives an optical depth of r 15 01 = 0.4 for the Fe XVII 15.01 A line. We adopt values 
of r 15 - 01 = 0.4 , fij = 2.66 , M = 55.8 amu, T = 10 6 - 8 K, n-Jn cl = n Fc xvu/n Fc = 0.63, 
n-ei/n-n = n Fc /n u = 1.6 x 10~ 5 , n H /n e = 0.8, n c = AA x 10 10 cm~ 3 . The iron abundance was 
obtained by matching the continuum of the predicted spectrum to the observed, while the 
electron density was estimated from the O VII triplet using the CHIANTI database. Using 
the above values we derive I ~ 8 x 10 8 cm = 8, 000 km. The uncertainty in the theoreti- 
cal line ratio is difficult to quantify but should be better understood once the discrepancy 
between the two laboratory measurements is understood. 
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Although the Fe XVII 15.26 A/15.01 A line ratio is a more reliable opacity diagnostic 
since both lines are populated by collisional excitation from the ground, and the laboratory 
data for this line ratio are in good agreement, (Laming et al. 2000), the 15.26 A line is 
extremely weak and falls on a bad pixel of RGS1. Testa et al. (2004) have found decrements 
in the O VIII Lya/Ly/3 ratio in the coronae of very active RS CVn binaries which are 
interpreted in terms of opacity. We were not able to do a similar measurement here as the 
Ly/3 line is blended in the RGS spectra. 



4. Conclusions 

While it is now accepted that resonance scattering can affect the flux of lines in the 
solar X-ray spectrum, especially when focusing on small active regions and flares, the results 
for quiescent stellar coronae remain controversial. Depending on the source geometry, line 
photons may scatter in or out of the line-of-sight and a net effect is difficult to detect 
when observing the integrated emission from the stellar surface. Our analysis of the X-ray 
spectrum of AB Dor has shown a significant increment in the Fe XVII 16.78 A/15.01 A 
line ratio during a flare event as compared to the quiescent state. Although the theoretical 
value of the 16.78 A/15.01 A ratio depends on the adopted atomic model for Fe XVII, the 
latest calculations that include both direct and indirect processes predict a decrease in the 
line ratio with increasing temperature (see Doron & Behar (2002) and Gu (2003)). This 
is opposite to what we observe in the flare. We therefore attribute this change to opacity 
effects. We demonstrate the potential of coronal opacity as a diagnostic and derive a path 
length of w 8,000 km (0.01 R*) for the coronal scattering layer. This value is two orders of 
magnitude greater than the path lengths determined in stellar transition regions and falls 
within the range found for solar coronal loops (Nakariakov et al. 2003). 

This work is supported by the Particle Physics and Astronomy Research Council (PPARC). 
M. Matranga would like to thank Queen's University Belfast for a research studentship. 
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Florence (Italy) and Cambridge (UK). We thank an anonymous referee for his/her comments 
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Table 1. Lines Used for the Emission Measure Reconstruction 





Ion 


logT max b 


Blends 


6.648 s 


Si XIII 


7.00 


Si XIII 6.688 


8.425 I 


Mg XII 


7.00 


Mg XII 8.419 


9.169 f 


Mg XI 


6.80 




10.238 i 


Ne X 


6.80 


Ne X 10.240 


11.171 q 


Fe XXIV 


7.30 


Fe XXIV 11.188 


11.426 q 


Fe XXIV 


7.30 


Fc XXII 11.442, Fe XXIII 11.441 


11.736 f 


Fe XXIII 


7.20 


Fe XXII 11.767, 11.795, Ni XX 11.832 


11.767 q 


Fe XXII 


7.10 


Fe XXII 11.795, Ni XX 11.832 


12.132 -f 


Ne X 


6.80 


Fe XXIII 12.175, Fc XXII 12.193, Ne X 12.137, Fc XXI 12.204, Fe XVII 12.123 


12.285 * 


Fe XXI 


7.00 


Fc XVII 12.264 


12.435 /9 


Ni XIX 


6.90 


Fe XXI 12.398 / , 12.422^, 12.490 9 , Fc XIX 13.462^, 13.464^, 13.521^, Ne IX 13.447-f 


13.447 • f « 


Ne IX 


6.60 


Fe XIX 13.425?, 13.462^, 13.464^, 13.521/, Fc XXI 13.507-^ 


13.534 f 


Fe XIX 


6.90 


tti v"T"v 1 n r r*ci tv vvt i o r t a tv vvtt -i o r oci tvt yv 1 o r r o 

be XIX 13.568, re XXI 13.574, le XXII 13.582, Ne IX 13.553 


13.699 9 


Ne IX 


6.60 


Fc XIX 13.735, 13.736, 13.795, Fc XX 13.736, Ni XIX 13.779 


15.176 " 


O VIII 


6.50 


Fc XIX 15.198, 15.221, O VIII 15.177 


16.007 /9 


O VIII 


6.50 


Fe XVII 16.006 /9 , Fe XVIII 16. 005-^, 16.072 9 , Fc XIX 16.027 9 , O VIII 16.006 9 


16.165 ' 


Fe XVIII 


6.90 


Fe XVII 16.240 


16.778 Sq 


Fe XVII 


6.70 




17.053 /9 


Fe XVII 


6.70 


Fc XVII 17.098 /9 


18.627 i 


O VII 


6.35 




18.967 


O VIII 


6.50 


O VIII 18.973 /9 


20.910 9 


N VII 


6.35 


N VII 20.911 


21.602 / 


O VII 


6.30 




21.807 * 


O VII 


6.30 




22.101 * 


O VII 


6.30 


Ca XVII 22.114 


24.779 9 


N VII 


6.30 


N VII 24.785 


28.465 9 


C VI 


6.20 


C VI 28.466 


28.787 " 


N VI 


6.20 




30.427 9 


S XIV 


6.50 


S XIV 30.469 


33.550 q 


S XIV 


6.50 




33.734 


C VI 


6.20 


C VI 33.740-^ 



a -^prcd is the predicted wavelenght of the transition in A. 

b Tmax is the maximum temperature of formation of the line in K. 



Note. — In the "A prct j" column t and q indicate the lines we used for the EM reconstruction for the flare and quiescent 
observation respectively. In case where the line is selected for both observations, in the "Blends" column f and q indicate 
the blends present with that line for the flare and quiescent spectrum respectively. 
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Table 2. Observed Line Fluxes and Ratios for the Fe XVII 15.01 A and 16.78 A Lines. 

Theoretical Ratios are also listed 





Flux (16.78 A) Flux (15.01 A) 


Obs. Ratio 


Theor. Ratio 




[10~ 5 cnts cm~ 2 s _1 ] 






quiescent 


29.7 ± 1.5 49.1 ±1.8 


0.60 ± 0.04 


0.57W (600 eV) - 0.47^ 


flare 


39.1 ±2.6 54.8 ±2.8 


0.71 ±0.06 


0.54W (800 eV) - 0.47( 2 ) 



References. - 



W Doron & Behar (2002); ( 2 ) theoretical ratios from CHIANTI v4.0 
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PN lightcurve 200 sec binsize 




7.3540x10 7 7.3550x10 7 7.3560x10 7 7.3570x10 7 7.3580x10 7 7.3590x10 7 

Time [s] 



Fig. I. — The AB Dor light curve derived from EPIC pn using a bin size of 200 sec. During 
the observation two flares were observed. The flare and quiescent time intervals (FLARE 
and QUIESI+QUIES2 respectively) used to extract the RGS spectra are shown. Times in 
seconds since JD 2450814.0 . 
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Fig. 2. — Observed spectrum (solid line), predicted spectrum (dashed line) and continuum 
(dotted line) of the flare. The predicted spectrum has been calculated using the new atomic 
calculation of Doron & Behar (2002). The wavelengths grid of the predicted spectrum has 
been shifted by 0.1 A for clarity. The observed over predicted ratio for the Fe XVII 15.01 A 
line is lower compared to the quiescent state. 
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Fig. 3. — The quiescent spectrum of AB Dor. Symbols as in Fig 2 
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